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Abstract 

In this paper, the particles of quantum gases, that is, bosons and fermions are regarded as g-ons 
which obey fractional exclusion statistics. With this point of departure the thermostatistical relations 
concerning the Bose and Fermi systems are unified under the g-on formulation where a fractal approach 
is adopted. The fractal inspired entropy, the partition function, distribution function, the 
thermodynamics potential and the total number of g-ons have been found for a grand canonical g-on 
system. It is shown that from the g-on formulation; by a suitable choice of the parameters of the 
nonextensivity q, the parameter of the fractional exclusion statistics g, nonextensive Tsallis (q ^ l) as 
well as extensive (q=l) standard thermostatistical relations of the Bose and Fermi systems are 
recovered. 
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1. Introduction 

For thermodynamical systems which have high densities of particles at very low temperatures, the 
criteria of being classical n?i «1 where n is the density of the particles and /lis the thermal 
wavelength is no longer fullfilled. That means for very small temperatures as well as for high densities 
the interaction in real systems are no longer negligible and quantum effects are exhibited. The influence 
of quantum effects can be very well suited by the Bose and Fermi gas models where the corresponding 
well known Bose-Einstein (BE) and Fermi-Dirac (FD) statistics have fundamental importance for the 
understanding of the physical systems. Indeed, a great number of phenomena have been explained 
within the framework of BE and FD quantum statistics. Are they, however, unique and adequate to 
enlighten all the physical properties of the nature? In this connection, where a fractal space is involved, 
the quantum systems are expected to violate the standard properties and the nonextensive 
thermostatistics [1-23] is required. 

The second point to be raised is how one can reconcile the standard statistics with fractional statistics, 
e.g. para-statistics [24], intermediate statistics [25], fractional exclusion statistics[26-27], etc? Can BE 
and FD statistics be unified as well as to include many other statistics? 

Therefore, by adopting the fractional exclusion statistics the generalization for quantum gases has 
been done within the context of nonextensive thermostatistics the so-called Tsallis thermostatistics and 
finally the thermostatistics for bosons and fermions are unified. 

The proposed formalism basically relies upon two parameters, one is usually called the entropic 
index; q e R ,which takes care of the nonextensivity and it includes the standard extensive quantities as 
a special case when q is set to be unity and the other is the statistics parameter g which may have 
fractional values representing the ability of the occupancy of the single particle states by a g-on, with 
g=0 corresponding to usual bosons and g=l fermions [28-30]. 
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The outline of this paper is as follows; In sections 2 and 3, for the benefit of the readers the origin of 
the fractal inspired entropies and derivation of Tsallis entropy are stated. In section 4, the generalized 
entropy of g-on gas is obtained. In section 5, due to the great importance of the generalized mean 
occupation number for g-ons which includes the standard extensive BE and FD distribution as a special 
case where g is taken for bosons and 1 for fermions, the generalized g-on distribution is derived 
following two different methods. In section 6, the partition function, the thermodynamics potential, the 
total number of a g-on system within the nonextensive formalism are obtained. In section 7, the basic 
equations of the thermodymamic potential and the total number of the g-on system are calculated as an 
application, and their relation to the Bose and the Fermi system are exhibited. 

2. Measure, Spectrum of Fractal Dimensions and Fractal Inspired Entropies 

Consider a non-interacting quantum gas of N g-ons which is not in thermodynamical equilibrium. Let 
the system be in heat and particle baths, in other words, the total energy E and the total number of g- 
ons N fluctuate, but the total energy and the total number of g-ons are conserved on the average, that 
is, a grand canonical ensemble is established. The accessible states of the g-on system is given by the set 
of the occupation numbers \n l ,n 2 ,...n k ,...n K } where n k denotes the number of particles in state k [31]. K 
is the number of partitions of N particles. When a measure is done on the set 



K 



M = Tp q . 



(1) 



where P j is the probability of the i.th partition. 



Let us introduce the definition D q as the spectrum of the fractal dimensions 



D = Urn 

1 x^n 




(2) 



or it can be written in an exponential form as 



M =Q ~lim8 (q - l)D(q) . 



(3) 
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where 8 is the resolution related to the measure. From the relation between the entropy and the 
spectrum of fractal dimension, one of the fractal inspired entropies which is known as Tsallis entropy 
can also be obtained. The relation between the entropy and the spectrum of fractal dimension is given as 

S=-k R limln8 D (4) 

where k B is the Boltzmann constant. 

On the other hand, in view of Eq. (1) and Eq.(2),the spectrum of fractal dimension is obtained in 
terms of probabilities: 



d i lnft g = i ln£/* 
q q-l In 5 In 5 q-l 



(5) 



or 



1 i=i 



D a — (6) 

q q-l In 5 



since 



InM^lnQ^ £P«,--1 (7) 
Then when Eq.(6) is substituted into Eq.(4),Tsallis entropy is obtained: 

S q =~ - 1 (8) 
q-l 

where K can be regarded as the number of accessible states of the system.In the limit q — > 1 the well 
known Shannon entropy is recovered. 

S l =-k B LP i lnP i (9) 



3 



3.Nonextensive Entropy of a G-on gas 

The entropy of a quantum gas can be expressed in terms of the statistical weights. The relation 
between the statistical weight of a g-on gas which takes into account the statistics which it obeys and its 
thermodynamics is provided through statistical mechanics. From this point of view, entropy is a bridge 
between micro and macrophysics. Thus , the entropy of a g-on gas system has a very close relation 
with the statistics it obeys. 

In this section, the fractal inspired entropy of a non-interacting g-on gas system, which is not in 
equilibrium has been obtained in terms of the occupation numbers. 

The statistical weight of a g-on gas system which is not in equilibrium obeying the fractional 
exclusion statistics, is given by [27,28] 

k=i n k ! [g k -gn k -{l-g)\! 
There are exactly f2(g)ways to distribute the indistinguishable nk g-ons over the gk states in a certain 
cell, where gk is the number of states in group k, i.e. the degeneracy, nk is the number of particles in 
these states (occupation numbers), g is a statistics parameter < g < 1 , indicating the ability of a 
particle to occupy a single particle state, with g=0 corresponding to usual bosons and g=l to fermions 
[36,37]. 

In order to find the nonextensive entropy of a g-on gas fractal dimension D q is needed. By using Eq. 
(10) and following the generalizing technique in Ref[15]. 

°° uf-vf-n q 



D{g,q)=Y, J * * (ID 
( q -\ )lno 



£j (q-l)lnd 
is obtained for the g-on system, where 

u k =g k +{n k -\){\-g) (12 a) 

and 

v k =Sk ~8 n k -(!-<?)• (12 b) 



Eq.(l 1) can be rewritten as [15] 

D (g>l)=L8k "f~\ k ~" k (13) 
^ (q-l)lnd 



K uf-vf-n? 



where for simplicity 



"fc-> — ->l + njfc(l-g) ( 14a ) 
Sk 



v k ^^^\-gn k , n k ^^- (14b) 



Sk Sk 
are taken in place of the corresponding expressions. 
Using Eq.(12) , the fractal inspired entropy of g-on system is found: 

SM^, }*"*-'* (15) 

where Eq.(14 a), (14 b) are used. The moment order q, which appears in the entropy for the g-on 
system, thereafter will be called as entropic index; in case of q — > 1 limit Eq.(15) reads: 

K 

S (g,l) = k B 8 Jl + n k (l - g)]ln[l + n k (l - g)]- [l - gn k ]ln[l - gn k ]- n k lnn k . (16) 



Eq.(16) may lead to the usual entropy of the Fermi system S(l,l) and of the Bose system S(0,1) for g=l, 
and g=0 respectively. Thus 



S(g,l)=+k B ^g k [±n k lnn k ±(l + n k )ln(l + n k )] (17) 

k=l 

where upper signs should be taken for the Fermi systems and the lower signs for the Bose systems. 

4. Derivation of The Generalized Distribution Function of a G-on System 

We assume that the H-theorem is valid for the generalized g-on system whose entropy is 
S(g,q). Therefore, the H-theorem, can be used to find the distribution function of the g-on gas. In course 
of time, the entropy of the g-on system must incline towards a maximum in accordance with 
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Boltzmann's H-theorem, i.e. equilibrium state. The problem is to find n k such that the equation 

of entropy for the g-on gas has a maximum under the constraints [15] 

N q =tskK (18a) 

k 

E q ="LgkZ k n q k (18 b) 

which means that the total particle number and the total energy of the g-on gas respectively are 
conserved. The constraints are imposed for a macrocanonical ensemble where the g-on gas of the total 
particle number and the total energy may fluctuate by conserving these quantities on the average. This 
method of Undetermined Lagrange multipliers take the constraints (18 a) and (18 b) into account. 
For this aim, the expression 

F(g, q )=S(g,q)-aN q -PE q (19) 

is constructed. The values of a and (3 have been identified as f5 = ^ a = -fiji where T is temperature, 

[i is chemical potential and Boltzmann constant k B is set to unity.When Eqs.(15), (18 a) and (18 b) are 
substituted into Eq.(19),after a partial differentiation with respect to n k and equating 8F{g,q) to zero, 
the following equation is obtained 

1 1 1 - St + n k (1 - g)] q - 1 + g(l~ gn k r 1 - nf 1 - a(q - l>,f 1 - fi(q - l)s k nf }= . (20) 
k=l 

where we assume n k ,g k »l and Stirling formula is used. Since we have imposed two constraints by 
implementing two Lagrange multipliers, we may assume that here the variations 8n k are mutually 
independent. Then each coefficient in Eq.(20) must vanish: 

(1 -g)[l + n k (1 - gjf-i + g (l- gnk )<H _ _ a{q _ -p( q -l)e k nf 1 =0 (21) 
thus 

(1 - g)[l + n k (1 - s)]*" 1 + g(l- gn k = nf 1 + a(q - l>,f 1 +p{q-l)e k nf 1 (22) 
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Unfortunately, Eq.(22) does not allow us to find an analytical expression for n k . In order to have a 
solution for n k now let us write the left hand side of Eq.(22) in a more compact form, i.e. 

r=0 (23) 

[i+(i-2sKF 

where C^_j is the usual combinatorial term .The term in the parenthesis on the left hand side of Eq.(23) 
is approximated as : 

i(i-grM-i)v +1 ]~(i-2gr • (24) 

Then Eq.(22) reads 

[l + (l-2gKr 1 =nf 1 + a(^-l)nf 1 + ^(^-l) G , nf 1 . (25) 
Now the most probable distribution of the N q g-ons over the single states is solved: 

n k {g,q)= 1 j (26) 

[1-(1-^(g,-^)]" _1 +2^-1 

where a and /? are substituted into Eq.(26). n k (g,q) can now be interpreted as the number of g-ons per 

energy level which unifies the form of the generalized quantum distributions for bosons and fermions. 

It is obvious that, for g=l the generalized Fermi-Dirac n k (l,q) and for g=0 the generalized Bose- 

Einstein distribution n(0,g)are recovered [15]. It should be remarked that, n k (g,q) also includes the 



standard distribution functions for extensive systems that is n k (l,l), n k (0,l) n k 



,1 



v 2 J 



for quantum 



(fermi and bose) and classical systems, respectively. Eq.(26) may also lead to a generalized Planck 
distribution for u=0 and g=0 [36-37]. Standard Planck distribution is recovered for u=0, g=0 and q=l. 
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5. Generalized Grand Canonical Description of a G-on System 



Because of the great importance of the generalized mean occupation number of a g-on gas we want to 
rederive it using the generalized grand canonical partition function. Along this line, we introduce for the 
first time a generalized partition function Z(g,q) for grand canonical ensemble of a g-on gas within the 
factorization method which has beens developed by the authors of the present article[15]; 



Z(T,V,g,ll,q) = Y[ 
k=l 



where 



1 



n/c = "k = 2 



q)n k x k ] 



1 



(27) 



(28) 



is taken. 



Eq.(27) can be regarded as a unified partition function of quantum gases. It is obvious that for g=0 the 
generalized partition function of the Bose gas Z(T,V ,0,/J,,q) and for g=l the generalized partition 
function of the Fermi gas Z(T,V,l,/J-,q) are included in Eq.(27). Following the method in ref.[15] the 
generalized distribution function of a g-on gas can be obtained by using Eq.(27), then, one ends up with 
the same given by Eq.(26). 

6.Thermodynamics of the G-on Gas 

In order to obtain the thermodynamics formulations of the g-on gas <£>(T,V ,g,fl,q) is required .The 
generalized thermodynamics potential of the g-on gas,however , is obtained by incorporating the 
partition function with the relation 

®(T,V,g,n,q) = -k B TlnZ(T,V,g,z,q) (30) 
where summation is done on all of the quantum states. By taking into account Eq.(27) one derives 

1 



K 



&{T,V,g,pL,q) = -k B TY^ 

k=l 



1 + 



(Or 



(31) 



where 



®k =[i-(i-<7kf 1 + 



(32) 



In q — > 1 limit, Eq.(31) may be compared with Eq.(4.2) of ref. [27]: 



K 



^(T,V,g,n,l) = -k B T^ 

k=\ 



1 + 



e Xk + g -1 



(33) 



In the mean time, in view of Eq.(32), Eq.(26) can be rewritten in terms of (0 k as 



1 



(34) 



Since the generalized thermodynamics potential is 

&(T,V,g,n,q) = -pV 



(35) 



the pressure of a g-on gas is obtained: 



(36) 



) 



While the total number of particles is found as 

k k (0 k + . 



(37) 



The state of the gas is determined by Eq.(36) and Eq(37). The energy of the gas, however, is given by 



E q =Y, n k(s><ihk=Y,-^— 

k k^k+g 



(38) 



where Eq.(34) is taken into account. 

In order to calculate the thermodynamics potential, the total number of particles, internal energy, etc. 
of the g-on system, while including the correction due to the nonextensivity one may need the Taylor 
expansion of the distribution function of the g-on gas system. Therefore, the Taylor expansion of 
Eq.(26) is written: 



n 



(q,g)=n k (g,l) + ^(q-l)x k 2 n k 2 (g,l) + 0(q-l) 2 



(39) 



where 
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It should be note that all terms in Eq.(39) are written in the exponential forms which mathematically 
are favourable. Apart from the first term,the rest of the terms can be interpreted as a correction due to 
nonextensivity of the entropy, each of which lead to zero in the limit q — > 1 . 

As a first concrete application of the generalized fractional exclusion statistics we want to calculate 
the properties of an ideal gas of nonrelativistic indistinguishable g-on system. Our aim is now to 
calculate the generalized grand canonical partition function, or more simply, its logarithm 
Q(T,V,g,z,q ), then the grand canonical thermodynamics potential according to Eq.(30) is found. In 
view of Eq.(31), Q(T,V,g,z,q ) is given as 



Let one of the g-on energies e k are those of the energies of a free quantum mechanicle particle in a 

box of volume V. The chemical potential or fugacity z are not fixed, but the total number of g-ons are 
conserved on the average and z has to be determined from the Eq.(37). For a large volume the sum over 
all of the single-g-on states can be given as the single-g-on density states, where spin is not taken into 
account. 



Q(T,V,g,z,q) = lnZ(T,V,g,z,q)=Y,ln 1 + 



1 



(41) 



k 



(Ok 



D(g) = Ae 2 



(42) 



can be written, where 




h 



(43) 



Thus, the summation in Eq.(41) becomes 



CO 



1 




(44) 
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where Eq.(42) is taken into account. When the integration is performed by parts Q(l,V,g,z,q) is 



reduced to: 



Q{T,V,g,z,q)=\Apyin{g,q)d< 



(45) 



After the substitution of the Eq.(39) into Eq.(45) and then rearranging, 



Q(T,V,g,z,q) = -Ap 



e /2 . 1/ .02 



-rfe+-( ? -l)j3 2 z 



^2 c -/te 



o(l-z(l-2^f 



2^J 



^2 c -/te 



J o(l-z(l-2^)e-^f + ^ l(l-z(l-2g>-^) 2 



^2 c -/te 



(46) 



is obtained. 

The integrals are performed using the related formulae of p 326 and p 330 of ref[38] since the 
necessary conditions are met. Thus one concludes that, for a nonextensive g-on system 



Q{T,V,g,z,q) = ?-Ap\zF 5/ ( g ,z) + Uq-l)(3 2 z 
3 { 72 2 



F y {g, z)- 2^ / {g, z)+V 2 F v {g, z) 

72 72 72 



(47) 



where the definition 



k(l-2gff 



(48) 



Z (l-2g)p n £i k n ~ l 

is adopted in which T(n) is the usual gamma function. Then the state of the nonextensive system is 



determined by the following equation: 

pV = k B TQ{T,V,g,z,q). 



(49) 



It is obvious that the value of Q(T,V , g,z,q) for a bose system is obtainable by taking g=0 while for a 
fermi system by g=l. e = plays a special role in the bose system since it is not taken into account by 
the integrals, therefore one has to explicitly account for the terms in the sums in eq.(37) and eq.(41) in 
the transition from summation to integration. 
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In the thermodynamics limit, that is in the limit of infinite volume with particle density held fixed, 
particles in the ground state, having no kinetic energy do not contribute to the pressure, then ,in terms 
of the thermal wavelength, the pressure of a Bose system become, as it is expected, 

lrp OO k 

" = f£^ (50) 



where T 



r 5} -4K 



K 2 J 



3 is taken. 



2 2 

In a similar manner, the thermodynamical properties of a Fermi system follows immediately from the 
logarithm of the grandpartition function of a g-on system, which is given by Eq.(47). In terms of the 
thermal wavelength, the pressure of the Fermi system is 

k k=l fc/2 

where g s is due to spin. 

7.The Total Number of The G-on Gas 

For a large volume, in accordance with Eq.(37), the total number of g-ons, could be written as 

N(T,V,g,z,q)=Az] pe ^ de + Uq-l)Ap 2 z] ^ ~f ^ B \ 2 (52) 
{e^-zil-lg) 2 J o(l-z(l-2g>-^) 

assuming that the single number g-on states can be calculated in terms of integrals. The performation of 
integrals in Eq.(52) using the related integrals in ref [38], p. 326 and p 330 leads to: 



N(T,V,g,z,q)=A z F y2 (g,z)+^(q-l)Af5h\F 7/ (g, z )-2^F 5/ (g, z )+^ 2 F v (g, z ^ 



n /2 /2 



(53) 



where F n (g,z) is given by Eq.(48). This is the formula for the total number of a g-on system. If g=0 and 
g=l is substituted in Eq.(53), the formula for the total number of a nonextensive Bose system 
./V(T,V,0,z,g) and for a nonextensive Fermi system N(T,V ,l,z,q) are recovered respectively. 
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For the total number of bosons in the extensive bose system, however, the parameter of statistic g=0 
and the entropic index q =1 are taken in Eq.(53). Thus, for an extensive Bose system the total number 
of particles 7V(T, V ,0, z,l) including e = energy states is found to be : 

N{T,V,0, z,1) = AzF 3/ {0,z)+-^- (54) 

/2 l-Z 

Z 

where the last term N = represents the contribution of the energy level e = to the mean 

l-z 

particle number and F 3 / (0, z) is given by Eq.(48). 

72 

In a similar manner, for the total number of fermions of the extensive Fermi system the parameter 
of statistics g=l and the entropic index q =1 are substituted in Eq.(53). Thus, for an extensive Fermi 
system, the total number of fermions /V(T,V,l,z,l) is derived. 
The total number of fermions is find in terms of the thermal wavelength as 

N(T,V,l,z,l)=g s ^t }~ l \/ (55) 



*=i k 



where T 



v2y 



71 

is taken, and this is an expected standard result [36,37]. 



8. Conclusions 

In this study, a unified grand canonical description of the nonextensive thermostatistics of quantum 
systems have been achieved. This study has been initiated from the point of departure of bosons and 
fermions, which obey different statistics, could be regarded as a one-particle, simply, g-on which obey 
fractional exclusion statistics. Then using the fractional exclusion statistics where g-on number g is 
involved, and g might have fractional values as well, the thermostatistical relations which concerns 
quantum systems are unified, within the framework of nonextensive formalism which is introduced by 
C.Tsallis where entropic index q is involved. In this context for a grand canonical ensemble the fractal 
inspired entropy, the partition function, the distribution function, the thermodynamics potential, the total 
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number of particles, the state of a g-on system are obtained. It is clearly shown that for g=0 and for g=l 
the above mentioned quantities lead to corresponding quantities of a Bose and a Fermi system, 
respectively, while taking into account the nonextensivity of the systems involved. Furthermore, when 
the entropic index q is set to one, the corresponding quantities for Bose and Fermi extensive systems i.e. 
standard relations are recovered. 
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